Effective temperature of a dissipative driven mesoscopic system. 
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We study the nonequilibrium dynamics of a mesoscopic metallic ring threaded by a time-dependent 
magnetic field and coupled to an electronic reservoir. We analyze the relation between the (non- 
stationary) real-time Keldysh and retarded Green functions and we find that, in the linear response 
regime with weak heat transfer to the environment, an effective temperature accounts for the mod- 
ification of the equilibrium fluctuation-dissipation relation. We discuss possible extensions of this 
analysis. 
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Recent technical developments in nanoscience have re- 
newed the interest in quantum transport in mesoscopic 
systems 0. Heat transfer in these devices is a crucial 
issue which remains poorly understood. 

A recent theoretical result in the field of glassy dynam- 
ics is that an 'effective temperature', T^ff, controls the 
low-frequency linear response in the limit of small en- 
tropy production (long waiting-time, weak drive) [^|^. 
Te/j is defined as the parameter replacing the environ- 
mental temperature, T, in the fluctuation-dissipation re- 
lation (fdr) evaluated at low-frequencies 0, "6] and it 
has the properties of a temperature in the sense that it 
can be measured with a thermometer and controls heat 
fiows and partial equilibration |^ |^. A similar phe- 
nomenon was found in the relaxation of quantum spin- 
glasses and in the Coulomb glass 0. 

The question then arises as to whether the notion 
of Teff plays a role in understanding some features 
of transport in quantum mesoscopic systems. To ad- 
dress it we study the modification of the fdr in a 
very simple, exactly solvable, mesoscopic device that 
consists in a metallic ring connected by a lead to an 
external particle and thermal reservoir and driven out 
of equilibrium b y a threading time-dependent magnetic 
field [mil El El El- In this setting there is an elec- 
tronic dc current along the wire and a heat flow towards 
the environment. 

We model the metallic wire with a system of non- 
interacting spin-less electrons described by a one dimen- 
sional periodic tight-binding chain with length L — Na 
{N is the number of sites and a the lattice spacing) , hop- 
ping matrix element w — W/A and bandwidth W: 
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(1) 



The time-dependent phase (pt with ip = ^/{^qN) and 
$0 = hc/e accounts for the external magnetic flux that 
we choose to depend linearly in time, $j\/(t) = $i. The 
contact term between the lead and the ring is = 
—wia{c\ca+c\ci). We modcl the lead and reservoir with 



a semi-infinite tight-binding chain with hopping ampli- 
tude Wa = Wa/4:, bandwidth Wa, and spectral density 
Pa(w) = 4^1 -u^/W^ 9{Wo,-\uj\). We assume that the 
reservoir is in equilibrium at temperature T and chemi- 
cal potential fi and that its properties are not affected by 
the coupling to the small ring. The Hamiltonian of the 
full system is then H = H""a + + H°' El- We use 
a system of units such that h — ks = = ^- AH our 
results have been obtained using a ring with = 20 and 
W = 1, and a reservoir with /i = — 1 and — 4. 

The dynamics of this problem is amenable to an ex- 
act treatment within the real-time nonequilibrium for- 
malism. The retarded and Keldysh Green functions are 



(2) 



Gf^{t,t') ^ -^([c,(^),c+(^')]-)■ (3) 

The angular brackets indicate an average computed in 
the grand-canonical ensemble, i.e. using H — [iM in the 
statistical weight with M = J2iLi ^Ici- The exact evolu- 
tion equations read 15] 



- / dhGgit,h)^§Jh,t') = S,,Sit - t') , (4) 



^Gf^(t,t')-GUt,t')w,,{t') (5) 
dt, [Gf,{t,h)i:^^{h,t') + G^,{t,h)J:§^{h,t')] , 

with Wkj (f) = w e"^* ^k,j+i +'w e~*"^* ^k,j~i and the sum- 
mation convention over repeated indices. 

We call g^'^ the Green functions of the bath, 



g^ir) = -^e{T) l^pl^{u:)e 



(6) 



duj 



, _ tanh(/3c./2)ps;c(^)e— , (7) 



T = t - t', (3 ^ l/T and pg^cu) = pa{oJ + fi). These 
functions arc stationary since the bath is assumed to be 



in equilibrium. The effect of the lead and reservoir is then 
captured by a correction at the contact site i = 1 : 



(8) 



If one assumes that the magnetic field and coupling 
to the reservoir had been switched on in the far past 
in such a way that any transient behaviour dependent on 
the initial conditions died out at the working times {t, t'), 
Eq. ISJl reduces to 

Gf^{t,t') = J dhdh Gf,(i,ii)S^(ii,i2)Gi^(t2,i') 

that yields once G^ has been computed from 

Eq. Q [3. 

In equilibrium a model-independent fluctuation- 
dissipation relation (fdr) between the retarded and 
Keldysh Green functions holds. Calling Gfj^{t,uj) the 
Fourier transform with respect to r the fdr reads 



G^^it,Lo)^ tanh(/3c./2) {Gf'^{t,Lo) - Gf^{t,Lo)) (9) 

where, for later convenience, we kept explicit a depen- 
dence on the observation time t that does not exist in 
equilibrium. Out of equilibrium there is no reason why 
such a relation should hold, and the temperature and 
chemical potential of the system are not even defined. 
Our aim is to determine how the relation ^ should be 
modified to account for the Green functions of the ring. 
First, let us discuss the values of the parameters T, 



$ and Vl = 



\Wlc 



for which we may expect a sim- 



ple generalization of the FDR. To this end we summa- 
rize the dependence of the current on time and these 
parameters in Fig. [isl H^ . In the lower panel we 
show the time-dependence of the current Ji^i^i{t) = 
2eiRe[we-"f'* {c\{t)ci+i{t))] for ? = 1, $ 0.6, T 0.05 
and two values of the coupling to the lead, Vl = 0.9 (solid 
black line) and Vl = 0.1 (dashed blue line). The current 
has the period tb = 27r/<l> of Bloch oscillations. Their 
amplitude is diminished for large enough T and Vl and 
for small enough $ as well as for very large $(> W). 

Dissipation gives rise to a dc-component in the current, 
Jdc = 'Tb^ /o"^ dt Ji^i+i{t), that is related to the power 
at which energy is dissipated in the form of heat through 
P — Jdc^ as has been verified exactly 0|. In the upper- 
left panel in Fig. we show the dependence of Jdc on the 
strength of the coupling to the reservoir, Vl, for several 
values of the electromotive force (emf), and tempera- 
ture. At fixed temperature Jdc is a nonmonotonic func- 
tion of Vl with the boundary value s Jd c(O) — Jdc(po) = 
and Jdc Vl when Vl ^ 0, E| . For rather large 
values of the bias ($ = 0.6 open symbols) the maximum 
occurs at Vl ~ 0.5. For small values of the emf ($ — 0.1 
filled symbols) the first increasing part of the curve is 
squeezed towards very low values of Vl and only the sec- 
ond decreasing regime is visible. For fixed $, the higher 
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FIG. 1: (Color online.) Dependence of the current on time 
and various parameters. Upper-left panel: the dc current Jdc 
against the coupling to the lead Vl for $ = 0.6 (open blue 
symbols) and $ = 0.1 (filled magenta symbols). DiflFerent 
symbols correspond to T = 0.05 (circles), T = 1 (squares) 
and T = 10 (diamonds). The large circle and square indicate 
the values of Jdc that correspond to the data shown in Figs.|5] 
|2|and^ Upper-right panel: Jdc against $ for T = 0.05. Open 
symbols with dashed blue line and solid symbols with dashed 
black line correspond to Vl = 0.1 and Vl ~ 1, respectively. 
Lower panel: Ji2(t) for <E> = 0.6, T = 0.05, Vl =0.1 (solid 
black line) and Vl = 0.9 (dashed blue line). The period tb ~ 
10 is indicated with an arrow. 



T the lower Jdc, as can be checked by comparing the 
curves with different symbols. 

The upper-right panel in Fig. ^ shows that Jdc is in 
approximate linear relation with $ for small biases; de- 
viations appear at a value $,„ that increases with Vl, 
e.g. « 0.075 for Vl = 0.1 and w 0.2 for Vl = 1. 

We now turn to the detailed study of the Green func- 
tions that we parametrize as Gf^'^ {t, t — t). In Fig.|3we 
illustrate the non-trivial dependence on t and t = <— i' by 



tracing, with solid black lines, ReG^(r) for three values 
of the total time t equally spaced within tb ■ The two pan- 
els show the evolution on the intervals < r < 8 (upper) 
and 30 < r < 38 (lower) . Interestingly enough, all curves 
fall on top of each other for r < t*(^, (3, Vl) ~ 1.5, and 
they later deviate demonstrating the breakdown of sta- 
tionarity. This trend can be understood by noting that 
Gf\j [t, ljj) contains a time dependent structure within 
the spectral range of the free ring, i.e. in the interval 
1^ + A*! ^ W/2, while the high frequency part is domi- 
nated by the stationary spectral features of the reservoir 
which set the quick response in the time domain. 

Figure [21 also shows a qualitative study of the fdr. 
The dashed lines are the inverse Fourier transform of the 
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FIG. 2: (Color online.) Breakdown of stationarity and the 
FDR. Real part of Gii(t,r) (solid black lines) and the inverse 
Fourier transform of the rhs of Eq. (dashed magenta lines) 
as a function of r for three total times t equally spaced in the 
interval [0,ts]. l^z, = 1, $ = 0.1, T = 0.05 and J^c ~ 0.03 
(red square in Fig.0. 



rhs of Eq. for the same values of t. The companion 
curves do not match and the fdr does not hold. We 
note, however, that for r < t* « 1.5 the behaviour is not 
only stationary but the fdr holds as well. 

As Fig. 121 demonstrates, the total time dependence is 
very complicated. Instead of studying the modification 
of the FDR for each value of t we found it natural to work 
with the averaged Green functions: 



1 



(10) 



In Fig. O we show a quantitative test of the fdr in a 
'classical' regime (loosely) identified as the values of the 
environmental temperature T such that the classical limit 
of the fdr, 5^(r) = ~i/{2T)drgSiT), r > 0, holds for 
the bath. We compare the left and right hand sides of 
the inverse Fourier transform of Eq. iPJ averaged over t 
as defined in Eq. (|10(l . In both intervals shown the accord 
between the two curves is rather good, proving that the 
FDR approximately holds when the driving force is weak 
and the temperature of the environment is high in such 
a way that P is extremely low. 

In Fig. ^ we test the fdr in the 'quantum' regime, 
T = 0.05. The symbols represent a fit using the func- 
tional form in the rhs of Eq. Q with T replaced by 
Teff = 0.143. For short time-differences, say before the 
first minimum in the upper panel, the FDR holds - and 
the fit also falls on top of the original curves. Indeed, if 
/3e//w 1, tanh(/3e//w/2) w 1 and the Keldysh and re- 
tarded Green functions are in linear relation. For longer 
time differences the deviations are clear but the fit ac- 
counts well for the data. 




FIG. 3: (Color online.) Study of the fdr in the 'classical' 
regime, T = 10. Real part of {Gii){t) (dashed blue line) and 
the inverse Fourier transform of the rhs of Eq. ^ (solid black 
line) averaged over a period rs, see Eq. dU, as a function of 
r. The two curves are almost identical, apart from some little 
deviations seen, e.g., at r ~ 25.5. The symbols are obtained 
using the same functional form in the rhs of Eq. Q with /3 
replaced by a fitting parameter Teff — 10.02 ~ T. Vl = 1, 
$ = 0.1, Jdc ~ 0.002 (red circle in Fig.Cll. 



Within the accuracy of our numerical solution, the ef- 
fective temperature describes rather correctly the modi- 
fication of the FDR whenever the dc-current is in linear 
relation with the emf $ (see Fig.^. Out of the linear re- 
sponse regime one cannot match the two sides of Eq. 
by simply using a single- valued T^f / . This limitation par- 
allels the one observed in classical strongly driven glassy 
systems 0, 0] . This is reasonable since the notion of an 
effective temperature was proposed to apply to a regime 
of small entropy production only 0, Q and this, in our 
case, corresponds to low P = Jdc^ values. 
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TABLE I: Estimates of the effective temperature and dissi- 
pated power within the linear response regime. 

In all cases T^ff > T, as also found in glassy systems 
relaxing from a disordered initial configuration and in 
driven classical systems 0. An heuristic argumement 
that explains this result is the following. The equations 
that govern the dynamics of the Hamiltonian H are the 
same as those of an effective Hamiltonian where the elec- 
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FIG. 4: (Color online) FDR in the 'quantum' regime, T = 0.05. 
Re(Gii) (dashed blue line) and inverse Fourier transform of 
the rhs of Eq. Q (solid black line) averaged over a period tb 
against r. The symbols correspond to the latter functional 
form with a fitting parameter Te// ~ 0.143. $ — 0.1, Vl = 1 
and Jdc ~ 0.03 (red square in Fig.0. 

trons are coupled to modes with energy quanta $ 14]. It 
is natural to expect that such a coupling enables a mech- 
anism for the electrons to spread out in energy, which can 
be interpreted as an effective increment of the electronic 
temperature. More quantitatively, we found that T^jf 
increases with the power dissipated, P ~ Jdc^- Some 
rough estimates are shown in Table 1. 

Let us discuss some open questions. So far Tg// is 
a fitting parameter with a meaningful physical interpre- 
tation. It remains to be proven [2, |^ , whether we can 
associate it to a bonafide temperature for nonequilibrium 
quantum systems with fermionic statistics. 

We did not need to use a ^e/ / since /Lt basically controls 
the period of the oscillations of the inverse Fourier trans- 
form of the RHS of Eq. that coincides with the period 
of Re(G{^)(r), see Figs.Oland^ In other nonequilibrium 
settings an upgrading of /i to /ig/ / might be necessary. 

We defined T^ff via the Green functions that are not, 
however, directly accessible experimentally. A similar 
description should apply to the current fluctuations 
that can be measured in noise experiments |l8|| . 

Our problem bears some resemblance with a classical 
ratchet for which a Tg/ / could be defined in the strongly 
driven limit for long wave-length observables It would 
be interesting to check whether this holds for quantum 
mesoscopic devices too. Even if we have not found a site- 
dependent Tef f for weak forcings we cannot exclude such 
a dependence in a more general setting. 

Finally, it would be interesting to check the modifica- 
tions of the FDR using different time-dependent magnetic 
fields, including disorder in the hopping rates, consider- 
ing different dissipative mechanisms, etc. 
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